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We review known and discuss new signatures of high-intensity Compton scattering assuming a 
scenario where a high-power laser is brought into collision with an electron beam. At high intensities 
one expects to see a substantial red-shift of the usual kinematic Compton edge of the photon 
spectrum caused by the large, intensity dependent, effective mass of the electrons within the laser 
beam. Emission rates acquire their global maximum at this edge while neighbouring smaller peaks 
signal higher harmonics. In addition, we find that the notion of the centre-of-mass frame for a given 
harmonic becomes intensity dependent. Tuning the intensity then effectively amounts to changing 
the frame of reference, going continuously from inverse to ordinary Compton scattering with the 
centre-of-mass kinematics defining the transition point between the two. 
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I. INTRODUCTION 

The technological breakthrough of laser chirped-pulse amplification [H has led to unprecedented laser powers and 
intensities, the current records being about 1 Petawatt (PW) and 10 22 W/cm 2 , respectively. Within the next few years 
these are expected to be superseded by an increase of about an order of magnitude each, for instance at the upgraded 
Vulcan laser facility Q|. Up to three orders of magnitude may be gained at the planned 'Extreme Light Infrastructure' 
facility [3 L This progress calls for a reassessment of intensity effects in QED and the new prospects of measuring them 
(see e.g. J], Q and Q for discussions of strong-field physics at Vulcan and ELI). There is a plethora of strong-field 
QED processes, which may be roughly categorised into two classes; loop and tree-level processes. The former include 
strong-field vacuum polarisation, the real part of which describes vacuum birefringence Q (for a recent discussion 
see 0) while its imaginary part signals Breit- Wheeler pair production [9]. Summing all orders of these one- loop 
diagrams (in the low-energy limit) one obtains the Heisenberg-Euler effective Lagrangian llOl which in turn yields 
Schwinger's nonperturbative mechanism of spontaneous pair production from the vacuum [111 ]. The optical theorem 
and crossing symmetry relate these one-loop diagrams to tree-level processes such as perturbative pair production, 
pair annihilation and Compton scattering. 

It is well known that one-loop processes are of order h and thus of a genuine quantum nature, while tree level 
processes generically do have a classical limit. As a result, one can introduce two distinct parameters which characterise 
the different physics involved. The first parameter is the QED electrical field, 

2 3 

E c = —j- = 1.3 x 10 18 V/m , (1) 

first introduced by Sauter [l^] in his analysis of Klein's paradox [ljj]. The presence of Planck's constant, H, and the 
speed of light, c, show that E c originates from a relativistic quantum field theory. In an electric field of strength E c 
an electron acquires an electromagnetic energy equal to its rest mass mc 2 upon traversing a distance of a Compton 
wavelength, X e = h/mc. Hence, E c may be viewed as the critical field strength above which vacuum pair production 
becomes abundant. This is also borne out by Schwinger's pair creation probability given by the tunnelling factor 



p ~ cxp(— ttE c /E) ll| where E denotes the 'ambient' electric field one succeeds in achieving. Currently, this is 
E ~ 10 14 V/m implying a huge exponential suppression. The perturbative variant of the Schwinger process, i.e. 
the (strong-field) Breit- Wheeler process 

HQ 

was observed about a decade ago in the SLAC E-144 experiment 
15l |l6| . There a Compton backscattered photon pulse of about 30 GeV was brought into collision with the 50 GeV 
SLAC electron beam. The huge gamma factor (7 ~ 10 5 ) led to an effective electric field close to the critical one, 
E 1 = 'yE ~ E c , as seen by the electron in its rest frame. 

In this context a second parameter comes into play, the 'dimensionless laser amplitude', given as the ratio of the 
electromagnetic energy gained by an electron across a laser wavelength A to its rest mass, 

ao= — - (2) 

This is a purely classical ratio which exceeds unity once the electron's quiver motion in the laser beam has become 
relativistic. It may be generalised to an explicitly Lorentz and gauge invariant expression [l^ |. For our present 
purposes it is sufficient to adopt a useful rule-of-thumb formula expressing ao in terms of laser power [l8| , 

a 2 ~ 5 x 10 3 P/PW , (3) 

so that ao is of order 10 2 for a laser in the Petawatt class. SLAC E-144, on the other hand, had ao of order one, 
hence by modern standards was in the low-intensity, high-energy regime. As high energy implies huge gamma factors 
and fields close to E c this is also the genuine quantum regime. 

In this paper we will concentrate on the segment of the QED parameter space that has become accessible only 
recently, characterised by large intensities, ao 3> 1, and comparatively low energies, to <C mc 2 , typical for experiments 
with an all-optical setup. We will thus stay far below the Breit- Wheeler pair creation threshold and will have to 
consider a process that is not suppressed by either unfavourable powers or exponentials. A natural process that comes 
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to mind is a crossing image of the Breit- Wheeler one, namely strong-field Compton scattering where a high-intensity 
beam of laser photons 7^ collides with an electron beam emitting a photon 7. In this case one has to sum over all 
n-photon processes of the type 



+ njL —> e +7 



(4) 



The study of this process(es) has a history almost as long as that of the laser. Intensity effects were addressed as early 
as 1963/64 in at least three independent contributions, by Nikishov, Ritus and Narozhnyi 
Kibble 23] and Goldman 



19, 20 



21 



Brown and 



24| . These works are written from a particle physics perspective i.e. essentially by working 



15, 1 



out the relevant Feynman diagrams. For modern reviews of these development the reader is referred to 
Nikishov and Ritus in [22] pointed out that a 2 , is proportional to E 2 and hence the photon density n 7 . The precise 
relationship is 



he 2 



2 +3„ 



(5) 



where v = fku/mc 2 is the dimensionless laser frequency and ~X 3 n 7 is the number of photons in a laser wavelength 
cubed. As the probability for the process |J4| is proportional to a 2 ," ~ n™ it becomes nonlinear in photon density 
for n > 1 and hence is called nonlinear Compton scattering [22I ]. Somewhat in parallel, the same process has been 
considered by the laser and plasma physics communities, with an emphasis, however, on the very low-energy and 
hence classical aspects. The appropriate notion is therefore nonlinear Thomson scattering. These discussions were 
based on an analysis of the classical Lorentz-Maxwell equation of motion, typically using a noncovariant formulation 
and neglecting radiation damping. Some early references are papers by Sengupta [251 ] . Vachaspati [261 ] and Sarachik 
and Schappert 27J. Since then there has been an enormously large number of papers from this perspective, many of 
which are quoted in the concise review 28 1. 

The main intensity effect can indeed be understood classically, the reason being the huge photon numbers involved, 
X 3 n 7 ~ 10 18 in a laser wavelength cubed. Due to the quiver motion in a (circularly polarised) plane wave laser field 
the electron acquires a quasi 4-momentum given by 



2k -p 



k 



P + 1L 



(6) 



Hence, the electron acquires an additional, intensity dependent longitudinal momentum caused by the presence 
of the laser fields. It may be obtained as the proper time average of the solution p M (r) of the classical equation of 
motion with p^ — p^(0) being the initial electron 4-momentum and k^ — 



the lightlike 4- vector of the wave [2 



Historically, ^ was first found in the context of Volkov's solution 30( of the Dirac equation in a plane electromagnetic 
wave. Volkov explicitly wrote down the zero component q° while the generalisation ([6]) seems to be due to Sengupta 
(note added at the end of his paper (jlj : cf. also the textbook discussion in [32, Chapter 40]). Upon squaring q one 
infers as an immediate consequence the intensity dependent mass shift, 

to 2 -> m 2 = m 2 (l + al) . (7) 



Although first predicted by Sengupta in 1952 [il| (see also 23[ 29|) it has so far never been observed directly j^]. A 
central topic of this paper will be to (re)assess the prospects for measuring effects due to the mass shift 0. 

The paper is organised as follows. We begin in Sect, [ill by reviewing the coherent state model of laser fields, 
which provides the link between classical laser light and light quanta (photons) in quantum theory. We then describe 
scattering amplitudes between these coherent states in QED, and how they are generated by an effective action 
describing interactions with a classical background field. We illustrate this theory with nonlinear Compton scattering, 
in Sect. IIII1 and give a thorough discussion of the kinematics of the colliding particles. In Sect. IIVI we give a variety 
of predictions for both Lorentz invariant and lab-frame photon emission spectra. Our conclusions are presented in 
Sect. M 
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II. QED WITH CLASSICAL BACKGROUND FIELDS 

We first address the question which asymptotic in-state we should take to describe the laser field. In principle, 
we would simply take the multi-particle state containing the appropriate number of photons of laser frequency and 
momentum, encoded in the 4- vector k = (to, k). We are immediately faced with the problem of not knowing exactly 
how many photons are in the beam. Similarly, as we do not know how many photons will interact with, say, an 
electron during an experiment, we do not know what to take for the out-state. To overcome these problems we 
invoke the correspondence principle: due to the huge photon number in a high intensity beam it should be feasible 
to treat the laser classically, as some fixed background field. Formally, this is achieved by describing the laser beam, 
asymptotically, in terms of coherent states of radiation 34, 35, 36, 3?| • The coherent states have the usual exponential 
form 

_C*(fc)at(fc)|0>, (8) 

where is the photon creation operator, C ll (k) gives the (normalised) polarisation and momentum distribution of 
the photons in the beam and TV is the expectation value of the photon number operator (the average number of 
photons in the beam). As usual, the state is an eigenvector of the positive frequency part of since 



a IM (k)\C) = VNC IM (k)\C), (9) 

Expanding the exponential in ([8]) , we see that calculating S-matrix elements between states including coherent pieces 
is equivalent to a particular weighted sum over S-matrix elements of photon Fock states. Working with coherent 
states may also be thought of, physically, as neglecting depletion of the laser beam, i.e. taking the number of photons 
in the beam to remain constant 38, 3^]. There is a natural connection between classical fields and coherent states as 
these states are the 'most classical' available, having minimal uncertainty. The associated classical field is essentially, 
as we shall see, the Fourier transform of the distribution function C. To see this we turn to the calculation of S-matrix 
elements between coherent states. 

Consider some scattering process with an asymptotic in-state containing the coherent state C, and some collection 
of other particles. For reasons which will shortly become clear, we will summarise all those particles not in the 
coherent state by 'in', so that our state is | in; C ). Similarly, we take an out state of the form ( out; C | where we have, 
in accord with the assumption of no beam depletion, the same coherent state. In operator language, we are interested 
in calculating matrix elements ( out; C | § | in; C ) of the S-matrix operator, 



§ = Texp 



OO 

- % - J At Hj(t) 



(10) 



Here Hi{t) is the interaction Hamiltonian (in the interaction picture) and T denotes time-ordering. We now write 
the coherent state (JSJ as a translation of the vacuum state (see e.g. [40(), 

\C) = f c \0), (11) 
where the commutator of the translation operator and the photon annihilation operator is 

[& li (k),f c ]=C ll (k)fc. (12) 

Extracting the translation operator from the states 1 , we are left with ordinary asymptotic Fock states but with a 
modified S-matrix operator, 

(out;C*|S|in;C) = (out | f^ 1 §T C | in) . (13) 



1 Under the usual assumption of no forward scattering. For the photons, this requires C M (fc') = for any scattered photons of momentum 
k'. 
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From the definition (fTOjl of S, the effect of the translation operators is to shift any photon operator A^ appearing 
in the interaction Hamiltonian by (the Fourier transform of) C^{k) which we denote by A^(x). Hence, the fermions 
interact with the full quantum photon field and a classical background field, A^{x). 

To be precise, and switching to the more common LSZ language, S-matrix elements are given by the on-shcll 
Fourier transform of Feynman diagrams with amputated external legs, as usual, but where the Feynman diagrams 
are generated by the action 

S[A,A,ip,lp\ = J d*x - -F^F**" + [$ + ieA + ie4] - m) V . (14) 

This is almost the ordinary QED action, but the photon field in the interaction term is shifted by A^, explicitly given 
by 



A^x) = VA / d A k C^k) + c.c 



(15) 

fc 2 =0 



This potential gives the classical electromagnetic fields associated with the momentum distributions C^k). Note 
that only the interaction terms of the action are affected by the presence of the background field, following (|13|) . 
We therefore have a natural and quite elegant way to calculate - we do not need to directly add up the individual 
contributions of the infinite series of terms generated by expanding the asymptotic coherent state. Instead, we simply 
include a classical background in the action which contains all the information about the chosen asymptotic photon 
distributions. Following 2tJ, 41 1 these results can be summarised by 



(out;C*|S|in;C) = (out | f^ 1 ST C | in) = (out | S[A] | in) , (16) 

where, on the right hand side, the asymptotic states are ordinary particle number states, with no coherent pieces, and 
the photon fields in the S-matrix operator are translated by A^ . 

Briefly, the same result can be recovered entirely in the path integral, or functional, language, following e.g. 
42 . Chapter 9.2]. The construction of S-matrix elements between coherent states proceeds just as for elements 
between Fock states, but the asymptotic vacuum wavefunctional must be replaced by coherent state wavefunctionals. 
Ordinarily it is the vacuum which is responsible for introducing the ie prescription into the action and from there into 
the field propagators. A coherent state does this and more - it translates the photon field in the interaction terms by 
the classical field ([15]). recovering (TT6"|) . 

Note that the modified action (fTi|) remains quadratic in the fermion field. All effects of the background are therefore 
contained in a modification of the electron propagator. The result is that, in Feynman diagrams, the propagator 
becomes 'dressed' by the background field A^ which surrounds the electrons. The propagator will be represented 
by a heavy line as in Fig. [TJ and has a perturbative expansion in terms of a free electron propagator interacting an 
infinite number of times with A^, as represented by the dashed line. The Feynman rules of the theory are otherwise 



+ 



+ 



+ . . . 



FIG. 1: Perturbative expansion of the electron propagator in a background field. 



unchanged from QED - there is a single three field vertex which joins the photon propagator and two of the dressed 
fermion propagators. This background field approach is equivalent to adopting a Furry picture [431 ]. in which the 
'interaction' Hamiltonian describes the quantum interactions while the interaction with the background A^ is treated 
as part of the 'free' Hamiltonian. 

In general, the fermion propagator will have no closed form expression. Since an intense background will be 
characterised by numbers larger than one (such as the intensity parameter ao), a perturbative expansion in the 
background is not suitable. We can of course use a coupling expansion, but this leaves us with an infinite number 
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of Feynman diagrams to calculate for any process, even at tree level. Fortunately, for the backgrounds considered 
in this paper and discussed below, the electron propagator is known exactly, allowing us to treat the background 
field exactly. We will now illustrate these ideas by applying them to the process of interest in this paper; nonlinear 
Compton scattering. 



III. NONLINEAR COMPTON SCATTERING 



In this process an electron, incident upon a laser, scatters a photon out of the beam. Using the background field 
approach described above, we use the action (fl"4|) . which contains the effects of the laser, and take the asymptotic in- 
and out-states to be, respectively, 

\p,X), (p',X';k',e\. (17) 

The pair (p, A) give the momentum and spin state of the incoming electron, similarly (p' , A') describe the outgoing 
electron and (k' , e) are the momentum and polarisation tensor of the scattered photon. Only one Feynman dia- 
gram contributes to this process at tree level, shown in Fig. [2] Note that the analogous scattering amplitude with 
'naked' electrons, corresponding to spontaneous photon emission in vacuum, vanishes due to momentum conservation. 
Calculating the corresponding S-matrix element amounts to amputating the external legs and integrating over the 




FIG. 2: Nonlinear Compton scattering Feynman diagram using dressed electrons (subscript L). 



single vertex position. Amputating and Fourier transforming the electron propagator in a background field gives us 



the solutions of the Dirac equation in that background 19|, |23|, [24J, |29( . We will write these electron wavefunctions as 



^p\{x). The S-matrix element of the process in Fig. [2] therefore reduces to 



[p',X';k',e\S{A]\p,\) 



ik' ■ 



P '\ 



(x) 



V2\k 



= 1 ^ P x(x) 



(18) 



To proceed we need to pick a background field so that we can explicitly calculate the wavefunctions ^S p \(x) and 
therefore the S-matrix element (fT5)) . This is the focus of the next section. 



A. Plane waves and Volkov electrons 



We will model the laser by a plane wave, A^ = -4 p (fc • x), with k a lightlike four-vector characterising the laser 
beam direction. The electron wavefunctions in such a background, or 'Volkov electrons' 30(, are known exactly. The 
propagator is also known and may be derived either in field theory or using a first quantised (proper time) method 
ll| . For a textbook discussion see [32j, Chapter 40] . The Volkov electron is 



k ■ x 



* p a (x) =e-^ x exp I ^J— J d£ 2ep ■ A(0 - e 2 A 2 (0 



1 



2k -p 



(19) 



where p 2 — m 2 and u p is the usual electron spinor. 

To better understand this wavefunction we specialise from here on to the case of A^ being a circularly polarised 
plane wave of amplitude a, 



A^ = cos(fc • x) + a 1 ^ sin(fc ■ x) , 



(20) 
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where a 7 - • k = and a,- • 



-a 5jk- The electron wavef unction becomes 



fpA^) = exp 



iq ■ x 



ax -p . ,, . a 2 -p , v 
ze- sm(K • x) + le- cos(/e • x) 



k ■ p 



k ■ p 



(21) 



We have not given the explicit form of the spinor part; it is easily written down and not needed for the discussion in 
this section. The important effect is that the electron acquires the quasi 4-momentum q defined in ^ from the laser 
field with the intensity parameter ao given by 



,2 _ 



(22) 



Technically, the origin of the quasi-momentum lies in a separation of the exponent in (|19ll into a Fourier zero mode 
and oscillatory pieces, with the zero mode causing the momentum shift, p — > q. Inserting the wavefunctions (|2ip into 
(|18[) . and omitting the details of the calculation 21(, we find that the scattering amplitude is a periodic function with 
Fourier series 



(p',A';fc',e|SL4]|p,A) = 



1 



(2|fe'| 2E q , 2E q )V* ^ 



M(n) <5 (4) (q + nk-q' - k') 



(23) 



A discussion of the amplitudes M(n) may be found in 32|, Chapter 100]. We will give below the explicit form of 
the squared amplitudes summed over spins A, A' and polarisations e. We do not consider polarised scattering and 
ang ular distributions in this paper, though these topics are interesting in themselves and are discussed in, for instance, 



44 



The sum in (|23|) is not a coupling expansion, nor does it appear directly from an expansion of the coherent state into 
Fock states. Instead, the momentum-conserving delta function in the n th term implies that M(n) can be identified 
with the amplitude for an electron of momentum q and mass m* , absorbing n photons of momentum k and emitting 
one scattered photon of momentum k' , 

e.(g)+n7(fc)-*e*(g / )+7(A0, ( 24 ) 
as illustrated in Fig. [3J As pointed out in the introduction, these multi-photon processes are the origin of the name 




n7(fc M ) 



FIG. 3: The effective Feynman diagram describing the n th harmonic process; an electron of mass m* absorbs n laser photons 
of momentum ky, and emits a photon of momentum k' . 



'nonlinear' Compton scattering. It is simplest to use the language of quasi momenta to formulate the kinematics of 
([23]) as (|24l) is a process involving effective particles. The asymptotic particle kinematics may be reconstructed from 
the relation ^ between p and q. The processes with n > 1 correspond to higher harmonics. Note that the n = 1 
process is analogous to ordinary, 'linear' Compton scattering. It is possible to normalise such that one does indeed 
recover the Compton cross section at ao = 0. We will use this below as a reference cross section for experimental 
signals. 
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B. Kinematics — forward and back scattering 

We will now study the kinematics implied by the momentum conservation in (|23|) . finding an expression for the 
emitted photon frequency in terms of incoming particle data which generalises the standard Compton formula for the 
photon frequency shift. This will later be used when we predict the emitted photon spectrum. 

The delta function in (123p implies the momentum conservation equation 

q + nk = q' + k' , (25) 

where q is given by ([6]), q' being defined analogously with p replaced by p' . As k is light-like we have 

q-k=p-k , q' -k=p' -k . (26) 

It is useful to first discuss the kinematics in terms of the Mandelstam invariants 

s n = (l + nk) 2 = ml + 2nk ■ p> m» , (27) 
t n = (nk - k'f = -Ink ■ k' < , (28) 
u n = (nk — q') = ml — 2nk ■ p' . (29) 




FIG. 4: (color online) Mandelstam plot for nonlinear Compton scattering. Solid segments of dashed lines correspond to allowed 
u n and t n regions for each depicted value of s n . 



Recall that these are not independent as s n + t n + u n = 2ml. As each of them depends on the photon number 
n they will be different for each of the sub-processes (|2"4"|) . The physically allowed parameter ranges are displayed 
in the Mandelstam plot of Fig. |U For the n-photon sub-process, if s — s n is held fixed, the allowed t and u ranges 
(highlighted in red/full segments of dashed lines) are 

imin = "2ml - s n - rni/s n Umax = m-t/sn back scattering 

^max = "min = 2m^ — s n forward scattering 

Obviously, the allowed i-range increases with photon number n. 
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In order to find the generalisation of Compton's formula for the scattered photon frequency (thus abandoning 
manifest covariance) we square (|25p so that we may remove q' from the game via 

nk ■ q = k ■ q = q ■ k + nk ■ k , (31) 

where the second equality follows directly from (|25p . Using the definition (JB]) and (126|) we trade q for p, arriving at 
an equation in terms of the asymptotic, on-shell momenta 

(TTl \ 
n + On \k ■ k' , (32) 

2k -pj 

where k' 2 = and p 2 = m 2 . We will assume, in what follows, that the electron and laser meet in a head on collision. 
That is, incident momenta are 

fc" = w(l,n), p» = (E p ,-\p\n), |n| = l. (33) 

Primed (outgoing) quantities are defined analogously. For a head on collision the only angle in play is the standard 
scattering angle 8 of the photon, determined via n ■ n' = cos 8. The remaining scalar products become 

np = \p\ , n' ■ p = — \p\ cos 8 . (34) 

From now on we measure all energies in units of the (bare) electron mass, m. This introduces the dimensionless 
parameters 

v= — , 7 =^ = coshC, ^7=^=sinhC, (35) 
mm m 



where £ is the rapidity such that 



P = -fL = - 1/7 2 = tanhC . (36) 



E p 



Of course, (3 and 7 are the usual Lorentz factors characterising the frame of reference from the electron's point of 
view. f3 = 0, for instance, corresponds to the (asymptotic) electron rest frame. Using these definitions, equation ([3^]) 
may be rearranged to express the intensity dependent scattered photon frequency as 

Til/ 

V ' n ^ = l + K„(a o )e-C(l-cos0) ' (3?) 
Here, e -< » is the (inverse) Doppler shift factor for a head-on collision, 



e-<=7(l-/3) = ^f (38) 

Going back to (|57|) we see that all the intensity dependence resides in the coefficient 

K„(a ) =nv- /37 + ao7(l -0)/2 = nv- sinhC + a 2 e~ c /2 . (39) 

Standard ('linear') Compton scattering is reobtained by setting n — 1 and cto = (no intensity effects). In this case 
(f3"T)) and (|3^|) give back the ordinary Compton formula, 

(40) 



1 l + (i/-/37)7(l-/3)(l-cos(9) l + (i/-sinhC)e-£(l-cos0) ' 

So, technically speaking, the two intensity effects on the scattered frequency are the replacements (i) v — ► nv in the 
numerator and (ii) ki(0) — » K n (ao) in the denominator. Explicitly, the latter is 

v — 0-y — > nv - /?7 + ^7(1 - /3)/2 , or v — sinhC — > nv - sinh ( + a„ e _c /2 . (41) 
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The possibility of the incoming electron absorbing n > 1 laser photons may be interpreted, in a classical picture, as 
the generation of the n th harmonic, modulated by both relativistic and intensity effects. Using a linearly polarised 
beam the first few harmonics have indeed been observed experimentally by analysing the photon distribution as a 
function of azimuthal angle, ip. The second and third harmonics have clearly been identified from their quadrupole 



and sextupole radiation patterns 47 1 



For each harmonic number n, the allowed range of scattered photon frequencies v' n is finite. The boundary values 
of this interval (which is the i-interval in the Mandelstam plot Fig. [4|) correspond to forward and back scattering at 
6 = and tt respectively, 

nV 

v' n (0) = nv, u (tt) = . (42) 

l + 2K„(a )e <» 

The assignment of minimum and maximum depends on the sign of c„, 

K n > v'(tt) < v'(9) < nv red shift 'Compton' 

(43) 

K n < =>• nv < v' n (9) < v n (ir) blue shift 'Inverse Compton' 

So, if K n > 0, the allowed scattered photon energies v' n are red shifted relative to nv, the energy of the n absorbed 
laser photons. This clearly includes the case ciq = 0, 7 = 1 and n — 1 which describes Compton's original scattering 
experiment in the electron rest frame. In accelerator language, this case sees the laser fired onto a fixed electron 
target; the laser photon transfers energy to the target, so that the scattered photon is red-shifted (1/ < nv). 

On the other hand, if K n < 0, the scattered photon's energy is blue shifted from nv. The situation when the 
photon gains energy from the electrons is often referred to as 'inverse' Compton scattering. This is of relevance in 



astrophysics, for instance in the Sunyaev-Zeldovich effect 48|, |49|, |50( . A particularly simple and important scenario 



is provided by the backscattering of the laser pulses, 9 — tt, in the high energy limit (inverse Compton regime). We 
take 7 3> 1 so that ~ 27 and we assume K n < 0, whereupon the scattered frequency becomes, from {37]), 



<W = 9 z 7 - 5 ; > ( 44 ) 

n\ J 1 1 ~2 1 o„,, 1 1 „2 .1 A ' V ' 



3 2 ^ A^' z nv 
1 + al + 2nv e? _ 1 + ag + i^nv 

where the approximation is valid for high energy. In this regime one may distinguish between two different limits 



v' n {jr) = Aj 2 nv/al if A^nv < 1< (45) 
^nW =7 if l + ao< 4 7n^ (46) 

It is the former subcase which is typically realised 2 for optical photons (y ~ 10 -6 ) and moderate values of harmonic 
number n. Thus, as long as clq < 27, the back scattered frequency (|45[) is (i) blue-shifted with respect to the incoming 
n th harmonic frequency nv and (ii) for n = 1, red-shifted compared to the linear 'kinematic edge' (the maximal, 



back-scattered frequency, v' max ) as emphasised already by McDonald [18|. Explicitly, this red-shift is 



Aj 2 v — ► ^ vla^ 1 1 > 1 . 4 7n^ « 1 « a„ . (47) 

From the definition of K n given in (|39|) it is clear that, given any fixed experimental setup (i.e. incoming electron 
energy and intensity parameters £ and do), K n will eventually become positive, and remain so for all higher harmonics 
with 

sinhC — ape^^/2 



n > 



no , (48) 



where denotes the nearest integer less than or equal to b. Thus, for a given experimental setup, scattered photons 
corresponding to harmonic generation with n > no can only have energies red-shifted relative to the energy nv 



2 SLAC E-144 had ■yu = O(l) so all terms in the denominator of 1144 l l were of comparable magnitude. 
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absorbed by the electron. Alternatively, we can fix n and so define a critical intensity, from the vanishing of K n , which 
allows us to tailor the emission spectrum. The critical intensity parameter is 

a o,critW = 27(1 + /3)C#7 - nv) = 2e c (sinhC - nv) = e 2C - 2nv e c - 1 > . (49) 

For ao = ao,crit(no) all harmonics with n > no (n < uq) will be red-shifted (blue-shifted). For the extreme choice 
of no = 1, all scattered frequencies will be red-shifted for intensities above ao jC rit(l), as in, for example, fixed target 
mode (7 = 1). We are, however, more interested in the colliding mode (high energy). Then, for 7 > 1, we can 
approximate a,g crit from (|49|) as 

<crit - 47 2 ~ 4fyrw . (50) 

When ijnv < 1 as above, ao,crit becomes effectively n-independent 

ao,crit ~ e c ~ 27 . (51) 

As a numerical example consider the facility at the Forschungszentrum Dresden-Rossendorf (FZD) with a 100 TW 
laser and a 40 MeV linac [5l|. This implies 7 = 80, c = 2x 10~ 6 and ao ~ 20, so that all harmonics are relatively 
blue shifted up to n ~ 3.9 x 10 7 - as we will see, emission rates at this n are basically zero. In this case, the critical 
value of ao, above which all harmonics (n > 1) are relatively red shifted compared to nv, is ao > 2"f = 160, an order 
of magnitude above the expected available intensity. One may verify, for example, that for ao = 200, K n > for all n. 

The discussion above will be illustrated in the next section when we discuss the photon spectra as a function of 
scattered frequency, v . In particular, we will see that, even if backscattering does not necessarily maximise the 
scattered photon frequency, it nevertheless gives us the strongest signal for which to search experimentally, namely 
the red-shift of the Compton edge (parameters permitting). 

To better understand the different behaviours of the harmonics, it is useful to write n n in terms of lab frame 
variables. For a head-on collision (which we assume), say along the z-axis, all momenta involved are longitudinal. 
The total 3-momentum, call it P, is then given by 

P = nk + q = nk + p + q L = m (nv — sinh£ + Qq e~^/2) z = m n n z . (52) 



The lab-frame physics involved in a head-on collision (p = 
the three terms contained in n n , 

n\k\/m — 
\p\/m = 
\Qh\l m = 



— (Pj/v)k) depends crucially on the relative magnitude of 

nv , (53) 
sinhC, (54) 
a 2 e- c /2. (55) 



Consider again Compton's original experiment with an electron at rest and ao = 0. This corresponds to q L = p = 0, 
so the only 3-momentum is that of the single incoming photon which delivers part of its energy to the electron and 
hence is red-shifted. If we now increase the electron energy in the lab (using a standard or wake field acceleration 
scheme) this red-shift turns into a blue-shift (v' > v) as soon as \p\ > \k\ — mv. This happens exactly where the 
total momentum, P = k + p, changes direction from pointing in direction k to — k. Hence, at this particular point P 
passes through zero, which, of course, defines the centre-of-mass (CM) frame where there is no frequency shift at all, 
v' = v. 

If we now turn on intensity (ao > 0) the total momentum acquires an additional, laser induced, contribution q L 
along k. So, in fixed target mode large intensity will result in a significant enhancement of the Compton red-shift. 
If, on the other hand, we assume colliding mode with a blue-shift at ao = 0, then the q L contribution in P works 
against the 'influence' of p. As a result, the blue shift v' > v at zero intensity is reduced, resulting in a red-shift of 
the kinematical Compton edge (1/^). If ao is large enough this latter red-shift may completely cancel the inverse 
Compton blue-shift. Again, this happens when the total momentum P = k + p + q L vanishes («i = 0) i.e. in the 
'CM frame' which is now an intensity dependent notion as q L depends on ao- 
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If we finally allow for higher harmonics n > 1, with the total momentum becoming P = nk + p + q L , we can 
balance p by increasing a or n or both. The transition point, K n = 0, defines a 'CM frame' for the n th process. At 
this point, the range of the n th allowed harmonic collapses to a point, v' n {9) — nv, as the 9 dependence in (|5T|) drops 
out. Strictly speaking, this can only occur for at most one value of n, but neighbouring n's will still have rather small 
spectral ranges (see Fig. [S] below). 



IV. PHOTON EMISSION RATES 
A. Lorentz invariant characterisation 

The S-matrix element represented by the Feynman diagram of Fig. [21 and given implicitly in (f23l) may readily be 
translated into an emission rate 1^> 32 1 . The non-trivial contribution to the differential rate for emitting a photon 
of frequency u>' = mv' per unit volume per unit time, in the n th harmonic process, i.e. the process (|24[) . comes from 
the differential probability 3 [il] 

^^jihf Uz{x)h n - l > (56) 

where x is the dimensionless Lorentz invariant 

k * k t j i 



> . (57) 



k ■ p' u n - 

The kincmatically allowed range for n th harmonic generation is given by the interval 

< x < y n , (58) 
_2nk-p s n 

Vn = k— = — T - 1 > , (59) 

which corresponds to the i-range given in (|30p. highlighted in Fig. |U The endpoints x — y n are located on the 
hyperbola su = rrQ. For x outside of this range the n th partial rate vanishes. 
The function Zn is 



4 , x 



,•2 



Z n {z) = ~Jl{z)+\2+^—\ jf l _ l {z) + Jl +l {z)~2Jl{z) , (60) 



the J n being Bessel functions of the first kind. Their argument is another Lorentz invariant 



2a x(y n - x) 

z(x) = J 2 . (61) 

yi V 1 + a o 



Both upper and lower limits of x correspond to z = and hence zeros of 3„(z) for all n > 1. The first few partial 
emission rates for E p = 50 MeV, u> = 1 eV (hence 7 = 10 2 , i/ = 2x 10 -6 ) and ao = 20 are plotted in Fig. [5] Linear 
Compton (ao = and n = 1) data is presented for comparison. 

The figure clearly shows the appearance of higher harmonics [n > 1) with, however, a reduced signal strength as 
compared to the fundamental frequency. Writing the Compton edge (|59"j) as 

77, 

Vn = y n (a ) = yi(0) 2 , (62) 

1 T On 



3 We normalise such that we recover the Klein-Nishina cross section for linear Compton scattering for n = 1 as ao — ► 0, see e.g. [32| . 
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X 

FIG. 5: (Color online) Partial emission rates (n — 1 ... 4) for nonlinear Compton scattering as a function of the Lorentz 
invariant x at intensity ao = 20, compared to linear Compton scattering (ao = curve). Horizontal log scale. 



where yi(0) corresponds to linear Compton scattering, we see that the edge x = yi(ao) of the first harmonic will 
always be shifted to the left by a factor 1/(1 + Oq). The same is true for the higher harmonics until n > 1 + a^. For 
ao 3> 1 these large harmonics will, however, be invisible due to their very small signal strength. 
To obtain the total rate one just sums over photon numbers n, i.e. over all harmonics, 

aw = ~ dw^ 

da; dx ' 

n— 1 

where it is understood that the n-th term is supported on < x < y n , with x given in (|57[) . The partial sums up to 
n = 30, 60 and 100 are shown in Fig. [HI along with the linear Compton spectrum. Again we note the significant shift of 
the fundamental Compton edge at x = yi(ao) together with side maxima due to the higher harmonics. Interestingly, 
the fundamental [n — 1) signal gets amplified due to superposition of the higher harmonic rates from Fig. [5j This 
suggests that, for ao > 1, the signal to noise ratio may become larger than for the linear case, while the full width 
at half maximum may become smaller. By tuning ao to an optimal value one may thus design X-rays of a given 
frequency and width. 



B. Lab kinematics: energy dependence 

Any actual Compton scattering experiment will be performed in a lab (frame) with the electrons either at rest (fixed 
target mode) or in motion. In what follows, we will assume the latter together with a head-on collision between laser 
pulse and electron beam (collider mode) as discussed in the previous section. In this case the kinematic invariants x 
and y n from (|5T[) and (|59|) become functions of the scattered frequency v' and the scattering angle 6, 

{1 - C ° Sdy (64) 



ef- (1- cos6V ' 
2n v 



(65) 
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FIG. 6: (Color online) Sum of partial emission rates from n = 1 . . . 30 (red, lower curve), 60 (cyan, middle curve) and 100 
(black, top curve) for nonlinear Compton scattering (head-on collision) at intensity ao = 20. The curves are indistinguishable 
for x ~ 10 -5 . Linear Compton data (blue, n = 1, ao = 0) added for comparison. 



Either the scattering angle 9 or the frequency v' n may be eliminated via (J37J) , allowing us to plot the emission rate as 
a function of v' or 9 respectively 4 . In this subsection we focus on the v' dependence of the partial and total emission 
rates which are depicted in Fig.s[7]and[8l respectively. Similar plots (for ao of order one) have been obtained before 



16, 



18, 



4J, [46( . Analytically the partial rates are 




FIG. 7: Individual harmonic spectra (n = 1, . . . 4) for nonlinear Compton scattering at intensity ao = 20 compared to linear 
Compton scattering (n = 1, ao = 0), as a function of v' . 



4 The relationship between angle and frequency spectrum l)37|l is invertible "provided K n ^ 0. For K n = the n th harmonic spectral range 
shrinks to a point (see below). 
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FIG. 8: Theoretical photon spectrum (sum of first 50 harmonics) for nonlinear Compton scattering at intensity ao = 20 
compared to linear Compton scattering (n — 1, ao = 0), as a function of v' . 



dW n = dW n dx = 3n{z) , gg . 
dv' dx dv' n n 

The allowed range for v' is given in (|42p and (|43|) . The argument z defined in (|6"Tj) becomes a function of i/' via its 
dependence on 

. nv — v' , . 

x = x n (v') = — , 67 

k 71 — nv + v' 

upon eliminating 6 from (|64[) via (|3T[) . 

For the parameters chosen (7 = 10 2 , v = 2 x 10~ 6 and ao = 20) Fig. s [7J and [5] are fairly similar to their invariant 
pendants, Fig. s [5] and [6] In particular, the previous shift in x now corresponds to a red-shift of the linear Compton 
edge by a factor of 1 + a^ ~ 400 from about 40 keV to 0.1 keV, i.e. from the hard to the soft X-ray regime. Note 
that the frequency range is still blue shifted relative to the incoming frequency v (corresponding to the left-hand 
edge in Fig. s [7] and [8] given by v = 2 X 10~ 6 ). Again, there is a noticeable enhancement of the total emission rate 
at v' n {ir) ~ A^v/a^, cf. (|45|) . due to the generation of peaks corresponding to higher harmonics, n > 1, with the 
peak values decreasing rapidly with n. We note that the edge values of the higher harmonics which are clearly visible 
in Fig. [7J get washed out by the superposition of more and more partial rates dW n in Fig. [5] This will reduce the 
visibility of the associated maxima, as will, of course, all sorts of background effects which have not been included in 
the theoretical analysis above. 

The properties of the photon spectrum depend crucially on electron parameters ((3, 7 or £) characterising the 
lab frame and, in particular, the intensity parameter ao. To illustrate this dependence along with the discussion of 
Subsection IIII Bl we have calculated the photon spectra as a function of ao, ranging from ao = 20 up to 300. The 
outcome is depicted in the movie-like sequence of plots of Fig. [9j As 7 = 100 the critical ao from (|49f defining the CM 
frame of the first harmonic is ao C rit(l) — 200 corresponding to the fourth plot in Fig. [21 There, the lower harmonic 
spectrum collapses to lines located at the individual harmonics with frequencies v' n = nv (marked by red vertical lines 
throughout). 

If we go through the whole sequence the following picture emerges. For small ao < docrit {n) , all harmonic ranges 
with counting label less than n are blue shifted. Plots 1 and 2 show the harmonic range for n = 1 (and part of 
n = 2), both to the right of their red end edges (y and 2v, respectively). The right-hand, blue end, maximum of the 
fundamental range is enhanced due to contributions of higher harmonics. For ao approaching its critical value the 
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FIG. 9: (Color online) Theoretical photon spectra for nonlinear Compton scattering for different values of ao (7 = 100) and 
incoming frequency u — 2 x 10 -6 . The vertical (red) lines correspond to frequencies nv. 

harmonic ranges shrink, and a gap between the first and second appears (Plot 3) so that the fundamental maxima 
become of equal height. At ao = aocrit(l) — 200 the first harmonic range shrinks (almost) to a point, with the 
neighbouring ranges also becoming very narrow (Plot 4). Once eto(l) becomes super-critical, all harmonic ranges are 
red-shifted (i.e. located to the left of the vertical (red) lines, v' n < nv), with the ranges increasing again and gaps 
closing (Plots 5 and 6). In Plot 6, the first and second harmonics overlap again, leading to maxima of different height, 
with the one at u[ = v being the larger. 

Thus, by tuning et we effectively change frames of references with <2ocrit(l) representing the border between inverse 
Compton scattering (blue-shift) and Compton scattering (red-shift). 
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C. Lab kinematics: angular dependence 



As mentioned earlier, the emission rates may be considered as functions of either scattered frequency u' n , or scattering 
angle 9, the two being related via (|37p . In terms of the scattering angle 9 the rates become 

dW n dW n dx e c xj 

In = ~dTdn = nu{i - C os9Y (i + Xn f Uzn) > 0<d <^ ( 68 ) 

where x n = x (for the n-th harmonic) and z n are to be viewed as functions of 9 (see below) . Our angular measure is 
df2 = dfsinf?, which is the solid angle measure up to a factor of 2tt, as the azimuthal angle <p does not contribute due 
to axial symmetry. Note that this is different for linear polarisation or, more generally, if there is another preferred 
direction which, for instance, could be induced by noncommutative geometry 52|. 

In terms of their angular dependence the various invariants may all be expressed, using (|37|) and (|57p . in terms of 
the variable x\ defined by 

2nf(l — cos#l 

x n {9) = n Xl {9) = —- K —— ^tt rr , 69) 

w w e<(l + cos(9) + e-C(l + ag)(l-cos6>) ' v ; 

with Xi between Xi(0) — and xi(ir) = yi as in (p)9")) , where 

<™> 

The argument of 3« in becomes 



z n (9) = n Zl (9) = 2n / U ° J r (l - r) , (71) 



where we have introduced the rescaled variable 



e~ c (l + a£)(l - cos9) „ 
^^ 1= eC (1 + co sg ) + e -C( 1 + ag)(l-cos^) ' (?2) 

As a result, zi becomes maximal for r = 1/2 and so zi is less than unity, 

which will be important later when we discuss the convergence of the emission rate sum. Solving r(#o) = 1/2 we find 
Z\ is maximised at the angle 

l+ag-e 2C 

6> = arccos 5 57 . (74) 

We will now relate these results to the emission spectra as functions of 9. In Fig. [TO] we show the angular distribution 
of the photon yield, as determined by (|68p . for the lowest individual harmonics, n — 1, ... ,5. For the parameters 
chosen (7 = 10 2 , v — 2 x 10 -6 and ao = 20) the largest signal is due to the fundamental harmonic, n — 1. This 
is also the only one contributing on axis, i.e. in the forward and backward directions, 9 = and it, respectively. 
For the classical intensity distribution this was also found by Sarachik and Schappert 27j. Thus, in particular, real 
backscattering at 9 = 7r only occurs for n = 1, while for the higher harmonics one has 'dead cones' with an opening 
angle of about 0.1 radians, slightly increasing with harmonic number n, as seen from the magnified plot in Fig. [TO1 
(right panel). 

The dead cones are controlled by the angle 9q from (f74|) : their opening angles are bounded by 9' Q = 9 — it. For 
1 « Bo < 7 2 the former are quite narrow such that most of the radiation (in particular the location of the maxima 
at #0) is near backward 5 . Quantitatively one finds that the dead cone opening angles are less than 

9' ~ a / 7 « 1 , (75) 



5 We mention in passing that the situation for linear polarisation is different. As pointed out by Esarey et al. [45l | for Thomson scattering 
with linearly polarised photons, odd harmonics do get backscattered (no dead cones). 
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FIG. 10: Theoretical photon spectrum for the first five individual harmonics as a function of scattering angle 9. Parameters: 
7 = 100, ao — 20; LEFT: vertical scale logarithmic; RIGHT: vertical scale linear, zoomed into range 2.7 < 8 < n. 



which, for the parameters of Fig. [10] corresponds to 0' ~ 0.2 radians. (For the intensity distribution of classical 
radiation the relation (|7"5"j) was found in [45jj.) 

To determine the total emission rate we have to sum (|68|) over all harmonic numbers, n. It is not entirely obvious 
that the ensuing series converges. To prove this we employ the Bessel function identity [s3 |. 



Jn±l(z) = ~Jn{z) T Jn(z) , 

the prime denoting the derivative with respect to the argument z, in order to rewrite Z in terms of J 2 and J^ 2 , 

2 / n 2 x? \ / 1 



3{z n ) = 2J„(nzi) 



1 + nx\ 



1 



+ 2J' n \n Zl ) 2 + 



1 + nx\ 



(76) 



(77) 



According to (|6"8")l . in the rates this is multiplied with an n-dependent factor n/(l + nxi) 2 . Thus, upon summation, 
we encounter series of the form 



n>0 v 



,N 



J 2 (nzi ) and 



^ (l + nxi) M 



where JVe {1,3} and M € {2, 3}. We can easily bound these series from above, for example 



£.71 



n>0 



(1 + nx\)'- 



■ J^inzi) < 22 nJn(nzi) = 5i , 



n>0 



(1 + nxi)'- 



n>0 

' Jn( nz l) < ^2 n3j n( nz l) = S 3 
n>0 



(78) 

(79) 
(80) 



(and likewise for J' 2 )- The series S\ and S3 on the right hand side are examples of Kapteyn series 54[ which are known 
to converge. Remarkably, some also have analytic expressions for the sum. These results do not seem particularly 
common, so we collect them in an appendix. Although we have not yet been able to explicitly perform our sums 
(which have a more complicated n-dependence than the Kapteyn series) we can now be confident that they converge. 
This is an extremely satisfying result confirming the validity of the background field picture we have employed and 
our analysis based around the summation of individual harmonics. 

Lerche and Tautz in [55| state that a summation of the first 1000 terms in Kapteyn series like (|79p or ([50)1 yields 
errors below 10 -6 for z\ < 0.95. We need to include Z\ values closer to one where the convergence rate is at its lowest. 
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This occurs near the angle 8q defined in (f74|) . Increasing the maximum harmonic number from 5000 to 10000 yields 
basically identical plots except that the height of the narrow peak at 8 increases as shown in Fig. [11] (left panel). 
The maximum is indeed located at 8 = 8 = 2.94 (or 8' ~ do/7 — 0-2) as given in (|74[) and ([75]) . The shoulder near 
8 = tt (0' = 0) is entirely due to the fundamental harmonic (n = 1). 




FIG. 11: Theoretical photon spectrum. Parameters: 7 = 100, ao = 20. LEFT: vertical scale linear, harmonics summed up to 
71 = 5000 (full line) and n = 10000 (dashed line). RIGHT: vertical scale logarithmic, harmonics summed up to n = 5000. 



Finally, we again vary ao and plot a movie of the angular distribution for fixed 7 = 100 in Fig. 1121 The main features 
are (i) a propagation of the main peak from near backward direction (when ao <C 27) to near forward direction (when 
ao 3> 27) consistent with the formula (|74|) for 80, (h) the appearance of a double peak which (iii) becomes symmetric 
for a ~ 27 at an angle 8 = ir/2. The latter situation corresponds to cos8 = 0, hence 

ag = e 2C - 1 ~ e 2C ~ 4 7 2 , (1 < a 2 < 7 2 ) . (81) 

This latter value (approximately) coincides with the critical ao of (|51[) . The locations of the two peaks in the spectrum 
are plotted in Fig. [T31 along with the angle 8 given in ([74]) . as a function of ao- It is clear from this plot that the 
maximum value of Z\ corresponds to the local minimum between the two peaks. 



D. Thomson limit: emission rate and intensity 

At this point one should mention that thorough discussions of the intensity distributions employing classical radi- 



ation theory have appeared before [27J, |45j . It is useful to check that our quantum calculations based on the Feynman 
diagrams of Fig. [3] describing nonlinear Compton scattering reproduce the results for nonlinear Thomson scattering 
in the classical limit. According to Nikishov and Ritus the classical limit is given by 

2np ■ k , . 

Vn = — V- « 1 > ( 82 ) 

mi 

which is just the statement that is the dominant energy scale. Note that this can be achieved by having large 
ao and may be counterbalanced by large n. Hence, harmonics with sufficiently large harmonic number n will behave 
non-classically (if they are observable at all despite their suppression). As y n is the upper bound for x n l[52"j) may 
equivalently be formulated as 

x n < 1 , (83) 
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FIG. 12: Theoretical photon spectrum as a function of 9, harmonics summed up to n — 5000 for different values of ao (7 = 100); 
vertical scale linear. 



such that we may neglect x n = nx\ on the left hand sides of (|T§1) and ([50)1 which hence coincide with S\ and S3 
in the classical limit. Even if (|83p no longer holds (i.e. for large n) contributions to the sum are still suppressed 
by J\. Comparing the quantum and classical (Compton vs. Thomson) rates by evaluating all sums numerically the 
graphs are indistinguishable. Plotting the relative difference for our parameter values one finds a small discrepancy 
near 9 = do, of the order of 1% (see Fig. Q3J). Note that the classical series S\ and S3 have a slightly slower rate of 
convergence (in particular near z\ = 1, i.e. 6 — 6q) where the suppression is mainly provided by J„(nzi), hence least 
efficient at z\ = 1. We have found for instance, that the peak in Fig. [14] increases from 0.4% to 0.7% when we increase 
the maximum n from 5000 to 10000. Nevertheless, Fig. [T4l provides a nice confirmation that for high intensity optical 
lasers the background can indeed be treated as classical to a very good approximation. 

We are left with relating photon production probabilities dW n to intensities d/„. This problem has also been 
addressed by Nikishov and Ritus [19( who state that the intensity is given by the zero component of the radiation 
4- vector, 



(84) 
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FIG. 13: The angular position of the maximum emission rates as a function of intensity, ao > 1 (blue/upper and lower curves), 
and the angle 9q which defines the maximum value of Z\ (red/middle curve). Harmonics summed up to n = 5000, 7 = 100. 




FIG. 14: Relative difference of photon emission rates |Compton — Thomson |/Compton as a function of scattering angle 9. 
Harmonics summed up to n = 10000, 7 = 100, ao = 20. 



We thus have dl n = mv'dW n or 



din 
dd 



= me 2 '' 



sin 2 6 a% (1 + nxi) 3 



Zn{nzi) 



(85) 



Compared to (|68[) we thus have an additional factor n/(l + nx\). In the classical limit, nx\ <C 1, this is just n so that 
([83)1 is bounded not by the Kapteyn series Si and S3, but by the analytically known series S2 and S4 as given in the 
appendix. 
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V. CONCLUSIONS 



In this paper we have (re)assessed the prospects for observing intensity effects in Compton scattering. The physical 
scenario assumed is the collision of a high-intensity laser beam with an electron beam of sufficiently high energy (7 > 
10 2 ) produced in a conventional accelerator or by a suitable laser plasma acceleration mechanism. In technical terms 
we were interested in the features present in cross sections or photon emission rates which are enhanced with increasing 
dimensionless laser amplitude, ciq = ea/m where a is the magnitude of the laser vector potential. The possible effects 
are of a mostly classical nature, being fundamentally due to the mass shift, to 2 — ► m 2 = to 2 (1 + a 2 ,) caused by the 
relativistic quiver motion of an electron in a laser field. Ranked in order of their relevance the main intensity effects 
are: (i) a red-shift of the kinematic Compton edge for the fundamental harmonic, u>' — 4j 2 ui — > A^oj/o^ for the 
parameters we have used, (ii) the appearance of higher harmonic peaks (n > 1) in the photon spectra and (iii) a 
possible transition from inverse Compton scattering (a/ > 10) to Compton scattering (us' < oj) upon tuning a®. The 
red-shift (i) may be explained in terms of the larger effective electron mass, to* > m, the generation of which costs 
energy that is missing when it comes to 'boosting' the photons to higher frequencies. This has, for instance, an impact 
on X-ray generation via Compton backscattering. To avoid significant energy losses (reducing the X-ray frequency) 
the amplitude ao should probably not exceed unity significantly. However, one is certainly dealing with a fine-tuning 
problem here, as item (ii), the generation of higher harmonics, improves the X-ray beam energy distribution. For 
ao > 1 there is a larger photon yield due to superposition of the harmonics and the full width at half maximum goes 
down. As a result, the X-rays tend to become more monochromatic once higher harmonics become involved. Item 
(iii), the transition from inverse to ordinary Compton scattering, once ao increases beyond 27 illustrates the energy 
'loss' just mentioned. When ao ~ 2j the lab frame can be interpreted as an intensity dependent centre-of-mass frame 
for which uj' n = ncu, at least for low harmonics. Thus there is no longer an energy gain of the emitted photons: the 
laser beam has become so 'stiff' that, in this frame, electrons begin to bounce back from it (gaining energy) rather 
than vice versa. 

The next step is to actually perform the experiments required for measuring the effects listed above. We emphasise 
that nonlinear Compton scattering provides a unique testing ground for strong-field QED as the process is not 
suppressed in terms of a or E/E c , by powers or exponentially. Hence, the experiments at Daresbury (7 ~ 50, ao — 2) 
[561 ] and the FZD (7 ~ 80, a ~ 20) planned for the near future should indeed be able to see the effects analysed 
in this paper. This will provide crucial evidence for the validity of the approach to strong-field QED adopted here, 
based on the electron mass shift, the Volkov solution and the Furry picture. 
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APPENDIX A: KAPTEYN SERIES 



The Kapteyn series 5J] (see also 57], Ch. XVII]) of the second kind involve squares of Bessel functions or their 



derivatives. We use the notation 

S N = E n N J 2 (n Zl ) , (Al) 

n>0 

S'n = E n N J' n 2 (n Zl ) , (A2) 

n>0 

(A3) 
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where < z% < 1 in keeping with our earlier discussion. The sums with a closed form expression are 

S- 2 = Y,n- 2 Jl{n Zl ) = Z l, (A4) 



4 

n>0 



SO^E J n( n ^ = 7==f - j > ( A5 ) 
n > V 1 

z 2 (4 + z 2 ) 



V- 4x2/ x z 2 (64 + 592z 2 + 472zf + 27zf) 



The first is a result of Nielsen [58| according to Schott who derived the second and third results [59|, p. 122], while the 
fourth can be found in [601 ] (note that our notation differs from that paper, which also contains a typographical error 
in their equation (24) for 5*2). The sums involving J' n are 



4 + 'iz\ 

1 u n [j izi) = — 

n>0 

4 T /2/ \ 64 + 624z 2 + 632z 4 + 45zf 
S' 4 = J] n 4 ^ (nzi) = 1 - ' 1 , (A9) 

n>0 ^ 1 ' 

given in j2?| and (r3o| . respectively. The latter paper also gives a double integral representation for the series S-i (there 
denoted F + ). Referring to a theorem by Watson [6l| the authors of 55j derive an iterative scheme for higher-order 
Kapteyn series, giving, for example, 

Si = (zi-^r) 5_i , (A10) 



1 — z{ \ dzi 

1 / _d_ 

1 — z\ \ dzi 



S 3 = - ^ ( zi— ) Si . (All) 
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